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Abstract 

In this paper we study the dynamics of the general case of Owen-Smith meta- 
physiological model, to explore the effects of seasonality on population fluctu- 
ations. The study will include the permanence, herbivore extinction, global 
asymptotic stability and existence of positive periodic solutions. Under certain 
assumptions, we obtained sufficient and necessary conditions which guarantee 
the permanence of herbivore and vegetation species, and existence of periodic 
solutions. The techniques we used here are, comparison method of differential 
equations, Lyapunov method and Brouwer's fixed-point Theorem. 
Keywords: metaphysiological, permanence, extinction, periodic solutions 

1. Introduction 

Owen-Smith in [1] presents the following metaphysiological model (a meta- 
physiological relates aggregated population dynamics to biomass gained from 
resources consumed, relative to physiological attrition and mortality losses [2, 
3, 4]), to explore how resource heterogeneity functions to dampen herbivore 
population fluctuations 
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where V, H represent vegetation, herbivore densities respectively, r is the max- 
imum relative growth rate of vegetation at time t, K is the vegetation carrying 
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capacity, i,„ is the maximum intake rate by unit herbivore, hi is the vegetation 
biomass at which herbivore intake rate reaches half of its maximum, bg is the 
half-saturation level which is the vegetation biomass at which the intake rate be- 
comes half of its maximum, w„ is the ungrazable amount of vegetation at which 
the intake rate of the herbivores becomes zero (vegetation reserve), C is the 
rate of conversion from consumed vegetation biomass into herbivore biomass, 
rup is the relative rate of physiological attrition in herbivore biomass, qs is the 
basic rate of loss of herbivore biomass through mortality due to senescence, and 
qo determine how the steepness of nutritionally-related mortality increases with 
diminishing food gains. For more background on system (1) see [1, 2, 5, 6]. 

Most environmental factors are highly variable and they show seasonal vari- 
ations, and in response, birth rates, death rates and other vital rates of popula- 
tions vary greatly in time [7]. In order to represent the effects of seasonality in 
his model, Owen-Smith [1] made the vegetation growth to be seasonal by using a 
step production function assuming that vegetation growth occurred continually 
during growing season, and then ceased during dormant season (i.e. r = r„ for 
the growing season weeks, and r = for the dormant season weeks). Seasonality 
enter Owen-Smith model (1) only via the maximum vegetation growth rate r. 
However food abundance and atmospheric temperature show high variation in 
between seasons, animal intake rates are responses to these factors. The vege- 
tation reserve represents the amount of vegetation that herbivores cannot eat, 
and this amount may also show distinct variations during wet and dry seasons. 
Herbivore death rate is caused by many environmental factors such as preda- 
tion, food abundance and food quality, and all these factors respond to seasonal 
changes. So, it is plausible to assume that all the parameters in the Owen-Smith 
model (1) change with time. 

To perform our analysis firstly, assume that all the model parameters are 
continuous w-periodic functions. Secondly, assume the half saturation rates for 
consumption and conversion are equal, i.e 6; = bg. Thirdly, and for biological 
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realism, we assume that w(io) > and /i(io) > 0. Finally, we simplify system 
(1) by making change of variables. Let H = h, V ~ v, a ~ ry, b ~ c = 

1, n, a = Ci,n, Pi bi~Vu, P2 bg-v^, 7 g^, p = w„, and R = m.p+qa+qs- 
With these assumptions the model given by (1) becomes the following: 

f = v{t) {a{t) - 6(t)z;(i)) - c{t)^§^^h{t) 

where a{t),b{t),c{t), f3{t),a{t),j{t), p{t), and are all continuous and posi- 
tive w— periodic hmctions. 

We developed the research method given in [7, 8, 9] and [10], and we estab- 
lished new sufficient conditions for which herbivore and vegetation having glob- 
ally asymptotically solutions and having periodic solutions. Also we established 
new necessary condition for which species for system (2) remain permanent. 

2. Preliminaries 

Before we going to our main results, we give some definitions, notations, 
lemmas and theorems that we are going to use in order to obtain the main 
results. 

Definition 2.1 

For a positive continuous and cj-pcriodic function /(i), we denoted by 
the average of f{t) over an interval of length lu 

^ Jo 

Definition 2.2 

For a bounded continuous function g{t) on Kwe use the following notations: 

= sup g{t), (?' = inf <?(<)• 
tem. ^^'^ 

The following theorem proves the positive invariant of any solutions to (2) with 

positive initial vales. Let ]R+ = [0,oo), and IR^j^ = {X eM^\X > O}. 
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Theorem 2.1 

If {v, h) e is a solution to (2) with positive initial values, then both 
will remains positive for all t > tg. 

Proof 2.1 

By integrating both side of system (2) with respect to t from to we obtain the 
following: 



v{t) — w(to)cxp 
h{t) = h{to)exp 



' \ . ,,,, . c{s) {v{s)-pis))h{s) 
a[s) — b[sjv[s) , ^ , ^ — I as 



vis) il3{s)+v{s)) 



a{s){v{s) -p(s)) 



-i?(s)-7(s)- 



{s) + v{s) 



/to V m+vis) 'v{s)-p{s)^ 
It is obvious when v{to) > 0, and h{to) > 0, both v{t) and are positive and 
then invariant for all t >to. The proof is complete. 

Definition 2.3 ([10]) 

The solution of system (2) is said to be ultimately bounded if there exists B > 
such that for every solution {v{t),h{t)) of system (2) there exist T > such 
that II {v{t), h{t)) II < B, for all t > to + T, where B is independent of particular 
solutions while T can depend on the solution. 

Definition 2.4 ([11]) 

The system (2) is said to be permanent if there exists a compact set F in the 
interior of IR^ = {(w, /i) g IF^Iw > 0, /i > O} such that all the solutions starting 
in the interior of IF^ ultimately enter F and remains in it. In other words the 
system is permanent if there exists positive constants S, and A with < (5 < A 
such that 



max 



lim infu(t), lim iid h{t) 

^ + 00 t— S- + 00 



lim supw(t), lim sup ft- (i) 



> s 

< A 



for all solutions of system (2) with positive initial values. 
Definition 2.5 ([12, 10, 13]) 

A bounded non- negative solution {v{t), h{t)) of (2) is said to be globally asymp- 
totically stable (or globally attractive) if, for any other solution {v{t),h{t)) of 
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(2) with positive initial values, the following condition holds: 

: 0. 



lim {\v{t)~v{t)\ + \h{t)-h{t) 

The following Lemma proves that, system (2) has a unique positive globally 
asymptotically oj-pcriodic solution (u*, 0). Using this theorem we can prove pos- 
itive invariance, permanence, and global asymptotic stability of the full model 
(2). 

Lemma 2.1 

If a{t) and h{t) are oj-periodic, and if Au:{b) > and Ai^{a) > 0, then system (2) 
has a unique positive uj-periodic solution (v* , 0) which is globaUy asymptoticaUy 
stable, where Ai^{-) defined by (3), and 



v*{t) = cxp<. / a{s)ds 



b{s) exp < / a{v)dv > ds + c 



(4) 



Proof 2.2 

When /i = 0, system (2) is reduced to the following differential equation: 

v{t) = v{t) [a{t) - b{t)v{t)] . 



(5) 



It is not difficult to show that (4) is w-periodic solution to (5). To prove v*{t) 
is asymptotically stable consider the Lyapunov function 



W{t) = \\n{v{t)} -\n{v*{t)}\ 



(6) 



where v{t) is any positive bounded solution to (5). Calculating the upper right 
derivative of W{t) along the solution of (5), we obtain: 



Wit) 
and then we have 



sgn{vit)-v*it)} 



v{t) v*{t) 



(7) 



W{t) = sgn{v{t)-v*{t)} [a{t)~b{t)v{t)-a{t) + b{t)v*{t))] 

= sgn{v{t)^v*{t)} [~b{t)v{t)+b{t)v*{t)] 

= -b{t) sgn{v{t)~v*{t)} [v{t)~v*{t)] 

= -b{t)\v{t)-v*{t)\ (8) 
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Assume that ( = max{| lii{w(t)}|, | ln{z;*(t)}|}. Let (* = exp {Q, then we have: 

\v{t)-v*{t)\ < C* W{t) 
substituting this back into (8) we obtain: 

W{t) < -b{t)C W{t) (9) 
integrating both side of (9) from to to t yields: 

W{t) <W{tQ)expl-C [ b{s)ds\. (10) 



And since (6) > 0, then hmt^^cx, W(t) = 0. Thereby Umt_>oo \v{t) -v*{t)\ = 
0. The proof is complete. 
Theorem 2.2 

If (w(t), h(t)) is any positive solution to (2). Tiicn there exists positive constants 
Ml and M2 such that 

lim sup{v(t)} < Ml, and lim sup{h{t)} < M2. (11) 

Proof 2.3 

From the first equation of system (2) wc have 

v{t) < v{t) [a{t) - b{t)v{t)] . 

Consider the following equation 

U{t) = U{t) [a{t) - b{t)U{t)] . (12) 

By Lemma (2.1), equation (12) has a positive oj-periodic solution U*{t) which 
is globally asymptotically stable. Let U{t) be any positive solution to (12) with 
U{Q) = v[Q). Using the comparison method of differential equations (see [14]) 
we have: 

v{t) < U{t), t > (13) 
since U*{t) is globally attractive, for given e > there exists a Tq > such that 

v{t) <U*{t) + e, t>To. (14) 
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Let Ml = maxt^[To,oo){U*{t) + e} and then 

lim sup{w(t)} < Ml. (15) 

To prove the second part of the theorem, from the second equation of system 
(2) wc have 

^^v m+'"*(t)+£ Mi~p{t)j ^ ^ 

where v*{t) is the periodic solution to (5). Integrating the both sides of (16) 
from To to t yields: 

hit) < hiTo) exp I J^^ ^-R{s) + j^—^^ j ds 

and since v*{t) is bounded over [Tq, t], then (-i?(i) + ^WI-'W+^-pW] _ ^^1^ 
is also bounded over the same interval. Now we have 

hit) < /i (To) exp {77} (17) 

where 

f' f , ^ a(g) [v*is)+e-pis)] ^is)Pis) \ 
Let AI2 = maxtg[2-g 00) {h (To) exp (77)}. Then we have 



lim sup{ft,(i)} < M2. 



The proof is complete. 



The vegetation biomass is always greater than its reserve Vu, and then we can 
conclude that limt_>oo hrf{w(t)} > «„. When Vu = 0, the following theorem 
prove that there exists an integer nii > such that Vnnt-^ao iiif{w(i)} > nii. 

Theorem 2.3 

Consider system (2). If (&) > and (a - <^^yP^ Af^^ > 0, where A^i-) 
defined by (3), there exists a positive constant mi < Mi such that: 

lim inf{v{t)} > mi. 
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Proof 2.4 

From Theorem (2.2), there always exists T2 > such that < mi < Mi. If 
v{0) > 0, from the first equation of system (2) wc have: 



v{t) > v{t) 
Consider the foUowing auxihary equation 
u{t) = u{t) 



, , C{t) (Ml - p(t)) _ , , ^ , ^' 

«W - ' , M2 - b{t)v{t) 



m 



m 



(18) 



(19) 



Let u{i) be the positive oj-pcriodic solution to (19), u{t) be any positive solu- 
tion to (19) with u{T2) = v{T2). From the comparison method of differential 
equations we have: 



v{t) > uit), 



t > To. 



From the global attractivity of u{t), there exists Ta > T2 such that 

t>T3. 



u(t) , , 



Let mi = min{^^}, t > T3, then 



lim inf{w(t)} > mi. 

f 00 



(20) 



The proof is complete. 

Theorem 2.4 

Consider system (2). If A^^ ( —R 



°'(^'~P) _ -7 ( 1 H I— 

13 + v* I \ ' V* - p 



> 0, where v*{t) 



is the unique positive u-solution to system (5), A^{-) defined by (3) and P{f) 
j3{t) + p{f), there exists a positive constant m2 < M2 such that 



lim iiii{h(t)} > 1112- 



(21) 



Proof 2.5 

We can choose eo > such that 



^<.($eo(0) >0, and A^{^,{t))>Q 



(22) 
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where e is same as in theorem (2.2) and 



1 



v*{t) -eo- pit) 

m 



(23) 



(24) 



From the first equation of system (2), consider the auxiliary equation with the 
parameter > : 



v{t) = v{t) 



(25) 



If A^{a) > 0, andA^ (b{t) + V^^(pj^ > Oj Lemma (2.1) tell us that system (25) 
has a unique positive w-periodic solution w,, (t). Let Vri{t) be any solution to (25) 
with Vri{0) — i'*(0), then for the above Eq a-nd from global asymptotic stability 
of w,j, there exists > T4 such that 



\v^{t) - Vrf{t)\ < ^, t> T5. 



(26) 



Let v*{t) be the unique positive w-periodic solution to (25) when t] = 0. From 
continuity of solution in the parameter 77, there exists rjo ~ ?yo(eo) and t (E 
[T5,T5 +a;] such that 



\v^{t)~v*{t)\ < ^, t&[T,,n+uj] 
Since w,,(i) and v*{t) are both w-periodic, then we have 



£0 



\v^{t)-v*{t)\<^, t>0, 0<77<7;o. 



Let < 771 < 770. From (28) we have 



(27) 



(28) 



(29) 



Suppose that (21) does not hold, then there exists Z e such that for all pos- 
itive solution {v{t, Z), h{t, Z)) of system (2) with initial condition (u(0), /i(0)) = 
Z , we have 

h{t, Z) < 27JI, (2r;i < eo), t > T5. (30) 
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On other hand theorem (2.2) tell us that, there exists Tq > T5 such that 

v{t,Z)<M,, Te > T5. (31) 
Applying (30) and (31) to the first equation of (2), it follows that for all t > Tq 



v{t,Z) > v{t,Z) 



a{t)-[b{t) + 2r^,^^v{t,Z) 



(32) 



Let uit) be a positive solution to (25) with rj — rji and u{Tq) — v(Tq, Z), then 
we obtain 

v{t,Z) > u{t), t > Tq (33) 

thus, u^i{t) is periodic solution to (25), so from its global asymptotic stability, 
for given e = there exists Tt > Tq such that 

v{t, Z) > u{t) > V.,, (t) - t>T^ (34) 

and hence from (29), and (34) we obtain 

v{t,Z)>v*{t)-e, t>T7. (35) 

From the second equation of (2) and by applying (35) we have 

h{t,Z) > h{t,Z) min{$eo(i),*eW} (36) 

integrating (36) with respect to t from T^ yields 

/i(i,Z)>/i(T7,Z)exp|^ min{$,„(s),$,(s)}ds|, t > T7. (37) 

Thus, from (22) it follows that h(t) — > +00, as t — > +00 which is a contradic- 
tion to the claim (30). So there must exists 7712 such that limt^oo inf h[t) > 7712. 
This complete the proof. 

3. Main Results 
Theorem 3.1 

If all the assumptions in theorem (2.2), (2.3) and (2.4) hold. Then system (2) 
is permanent. 
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The following theorem gives the necessary condition for permanency of the 
system (2). 

Theorem 3.2 

System (2) is permanent if and only if 

A. {-R + - ^) > (38) 

where Ai^(-) defined by (3), and v* (t) is tlic unique positive periodic solution of 
(12). 

Proof 3.1 

Suppose that 

A. (-i? + ^$9^ - ^) < (39) 

with this assumption we only need to show that limt i,oo h{t) = 0. From the 

first equation in (2) we have 

v{t) < v(t) [a(t) - b{t)v{t)] . (40) 

Consider the following auxiliary equation 

u{t) = u{t) [a{t) - b{t)u{t)] . (41) 

Let v*{t) be the unique positive w-periodic solution to (41) and u{t) be any 
positive solution to (41). From global asymptotic stability of v*{t) and using 
the comparison method, for any positive e < 1, there exists T'^^ > such that 

u{t) <v*{t)+e, t>T^^\ (42) 

By comparison method and from (42) we have 

v{t) <v*{t)+e, t>T^^K (43) 

It follows that 



1 r (- RU\ ^ U\ ^*it)+^-p{t) _ lit) iv*it)+e-m) \ 
u^Jo \ ^ ^ ^ m+v*it)+e v*{t)+e-p{t) ) 



dt < -eo- (44) 
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From (43) we also have 

- + «w ^^^^^ — wrwr) < ^''^ 

We will show that for the above e, there exists a T^^) > T^^) such that h(T^^^) < 
e, otherwise from (45) wc have 



< h{t) < /i(r(i)) 



exp 



T(i) V P(s) + V(S) 



(46) 



as t — ?► +c>o, which is a contradiction. 
Let 



, / , a(t)(i;*(t)+£-p(t)) jmm+s^m]\ 



(47) 



which is bounded for e e [0, 1]. Now we only need to show 

/i(r(3)) < eexp{M(£)a;} (48) 
otherwise there exists a T^'^) > T^^^ such that 

/i(t(3)) > £exp{M(£)cj}. (49) 

By the continuity of h{t) there exists T^'^^ e (r(2),T(3)) such that /i(T('*)) = e 
and /i(t) > e for t e (r(4),r(3)). 

Let P be the nonnegative integer such that Tt^) G (T^"*' + Plo, T*"*) + (P + l)w) , 
from (45) we have 

eexp{A/(e)w} < /i(T(3)) 

< MT^^^)) exp {/^: (-i?(.) + ^^Mlgh^ _ imgi)^) rft} 

- F exr. / r^'"+-^" + f 7?rn + "(*)('"(«)-p(«)) _ 7(t)[^(«)-/3(0] A j.l 

— fcexp|jj,(4, +Jt(-i)+Pw 1^ + 0(t)+i,(t) i'(t)-p(t) J j 

S £exp |jj,(4)+p„ ^^ n{T)+ i3{t)+v{t) v{t)-p{t) j"'^) 

< eexp{M{e)Lu} 

which is contradiction. It follows that h{t) — as t — y +oo, which contradict 
that h{t) is permanent, so (39) must be false i.e. (38) holds. This complete the 
proof. 
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The following theorem gives the sufficient conditions for (2) having a globahy 
asymptotically stable solutions. Let 



r = {{v,h) eW^\mi < V < Mi; m2 < h < Ah} . (50) 
Suppose that 

inf / h(f^ -Jr- p(f\ p(«)("H)-(™i^-p(«)<3(«))^ _ a(t)P(t) + a(t)p(t) 



t{t)m-i{t)p{t) 



> 0, 



{,n, - pit))iv{t) - Pit)) j ^ ' (-5^^ 

and 

■ r / m[-m^iP(t) + m)-i>it)pit)ifiit) + p(t)v{t))] \ „ 
mite]R|qrj M,«(t)(/3(t) + MjM,«(t)(/3(t) + «(*)) J ' 

where and A/^ defined by (50). 

Theorem 3.3 

Let {v{t),h{t)) be a bounded positive solution of (2). If (51) holds, and if 



Auj a 



(a - ^^^M,) > 0, [-R + - 5^) > 0, (52) 

where A(^(-) defined by (3). Then {v{t),h{t)) is globally asymptotically stable. 
Proof 3.2 

Let {v{t),h{t)) be any solution of (2) with {v{to) > O,h{to) > 0). The set F 
defined by (50) is an ultimately bounded region of (2), so there exists a T^^^^ > 
such that (u(t), h{t)) G F and {v{t), h{t)) e T for all t > to + T^^l Consider a 
Lyapunov function defined by 

X{t) = \\nv{t) -lnv(t)\ + \lnh{t) -\nh{t)\; t > to. (53) 



Calculating the right upper derivative of (53) produces 

v'it) _ (v(t))' 



h'jt) _ (hit))' 
h{t) hit) 



X'(t) =sgn{ln«(t) -ln{}(t)} - + sgn{ln ^t) - In ^t)} 

= sgn{ln.(i) - Invit)} [~bit)iv{t) - m) - c{t) { ^t^itT^W) " Iwffli 



-sg: 



n|ln/i(i) - In /i(t) I 



a(t)C»(t)-p(t)) _ 7(t)(/3(t) +»>(«)) _ a(t)(a(t)-p(t)) , 7(t)(;3 W + 
;3(t)+f(t) ''(t)-p(*) /3(t)+i)(t) «(*)-p(t) 
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< S! 



-^i^J v(t)vmm+v{t))m)+Ht)) ) - m\ j- 



< 



-b(t) - c(t) ^t*)^"'! + ~ - P(*)/3(t))^(*) I a(t)/3(t) + a(t)p(t) 

^ A/if)(t)(/J(t) + mi)(/?(i) +{)(*)) (/3(t)+mi)(/?(t)+f)(t)) 



7(t)/J(t)-7(t)p(t) 
("H -P(«))(«(*)-P(t)) 



- V{t)] - C(t) 



t.(t)(rni[/3(f) + f)(t)]--D(t)p(t)[/3(t)+p(t){)(t)]) 
M^v(t){l3{t) + M^){f3(t) + v{t)) 



h{t) - h{t) 



h(f) + r(f) + ~ ("'i^'^*) - Pit)0it))Ht) _ a(t)l3(t) + a{t)p(t) 

^' ^' Mi*(t)(/3(t) + mi)(/3W + «(t)) (/3(t) + mi)(/3W + *(*)) 



7(t)/3(t)-7(t)p(t) 



and since 



v(t)(7n^[l3(t)+v(t)]-v(t)p(t)[m + P{tXm 
M^v{t){p(t) + M^){l3{t) + ^(t)) 



h{t) - h{t) 



■ c I Uf\ , (f\ P(*)('"i + i'jt)) - {m^iit) - p(t)l3{t))h{t) _ a(t)fi(t) + a(t)p(t) 

teM\ \ )^A) Af^C(t)(0(t) + mJ(/3(t) + fi(t)) (/3(t) + mi)(/3(t) + C(t)) 



7(t)/3(t)-7(t)p(t) 



> 



(^l -p(t))(«(t)-p(*)) 

and inft ^ (c(i;) + ~ ^(')p(«)[^(*) + p(*)^(«)]) 



A/ii)(t)(/3(t) + A/i)(/3(t) + fi(t)) 

there exists a positive constant /U such that 



> 0, 



X'{t)<-^l(\v{t)-v{t)\ + \h{^)-h{t)\\, t>ta + T^'^\ (54) 
Integrating both sides of (54) with respect to t from T + tg yields 



X{t) - X{to + r(5)) <-fi (\v{s) ~ v{s)\ + his) - h{s) ) ds 



X{t) + ^l / {\v{s) ~ v{s)\ + h{s) - h{s) ) ds < X{tQ + T(^)) 



which yields 



M / f k(s) - «(s)| + /i(s) - Ks) ) rfs < x(tQ + r(5)) 

./to + T(=) ^ ^ 
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from here we can conclude that 



(\v{s) - v{s)\ + h{s) - h{s) )ds< n~^X{to + T^^') < +00 

and hence \v{t) ~ v{t) \ + h{t) - h{t) e ([to + T^^\+oo)). Since v{t),h{t) 
is bounded and vit), h[t) is ultimately bounded, this implies that they all have 
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bounded derivatives for t > ta + T . It follows that \v{t) — v{t) \ + h{t) — h{t) 
uniformly continuous on [<o + T^^\ +00) and hence using By Barbalat's Lemma 
[15] one can conclude that 



lim 

t< — hoo 



v{t)-v{t)\+ hit) -hit) 



0. 



(55) 



This shows that ^■O(t), is a globally asymptotically stable solution to (2). 
The proof is complete. 

Let (f (w, to, (uq, ho)) , hiuj, to, (wq, ^0))) be the solution of (2) through the point 
ito, ivo,ho)). 
Theorem 3.4 



Consider the set (50). If -R 



7/3 
Mi-p 



> 0, A 



and A^(a-^) 



(«(m1-p)) 



> 



> , where A^i-) is deGned by (3), and v is the periodic 
solution of (5), then system (2) has at least one positive td-periodic solution 
iv* , h*) which lies in T. 
Proof 3.3 

Define a Poincare shift operator a : — 5- , by 

cr iivo,ha)) ^ iviuj,to,ivo,ho)), hiuj,to,ivo,ho))) , (uo,^o)eIRt (56) 

Thus, from theorem (3.1) we know that F is positively invariant of (2), and hence 
(7(r) = r. Since the solution of (2) is continuous with respect to the initial value 
(to, (^Oj^o)), then the operator a is continuous. From the definition of the set 
r we know that it is bounded, closed, and convex. By Brouwer's fixed-point 
Theorem [16] we know that a has at least one fixed point 

iv*,h*) = aiiv*,h*)) = iviu,to,iv*,h*), /i(w,to,(«*,/i*)) (57) 

So there exist at least one periodic solution (u*, h*) G F. The proof is complete. 
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